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ABSTRACT 

We have observed 15 red giant stars in the relatively massive, metal-poor globular 
cluster NGC 4833 using the MIKE spectrograph at Magellan. We calculate stellar 
parameters for each star and perform a standard abundance analysis to derive abun¬ 
dances of 43 species of 39 elements, including 20 elements heavier than the iron group. 
We derive ([Fe/H]) = -2.25 ± 0.02 from Fe i lines and ([Fe/H]) = -2.19 ± 0.013 
from Fe ii lines. We confirm earlier results that found no internal metallicity spread in 
NGC 4833, and there are no significant star-to-star abundance dispersions among any 
elements in the iron group (19 ^ Z ^ 30). We recover the usual abundance variations 
among the light elements C, N, O, Na, Mg, Al, and possibly Si. The heavy-element dis¬ 
tribution reflects enrichment by r-process nucleosynthesis ([Eu/Fe] = -1-0.36 ± 0.03), 
as found in many other metal-poor globular clusters. We investigate small star-to-star 
variations found among the neutron-capture elements, and we conclude that these 
are probably not real variations. Upper limits on the Th abundance, loge (Th/Eu) < 
—0.47 ± 0.09, indicate that NGC 4833, like other globular clusters where Th has been 
studied, did not experience a so-called “actinide boost.” 

Key words: globular clusters: individual (NGC 4833) — nuclear reactions, nucle¬ 
osynthesis, abundances — stars: abundances — stars: Population II 


1 INTRODUCTION 


Globular clusters are the remnants of some of the most vio¬ 
lent epochs of star formation in the history of the Universe. 
The compositions of the present-day stars provide a detailed 
chemical inventory to probe those ancient nucleosynthesis 
and enrichment events. Nevertheless, detailed chemical in¬ 
ventories for some Galactic globular clusters are still rela¬ 
tively unknown. 

Only recently were the chemical compositions of a large 
sample of stars exam ined in one such cluster, NGG 4833 
llGarretta et al.l l2014l l. The properties and orbital param¬ 
eters of NGC 4833 are listed in Table [T] NGC 4833 has 
one of the more extended, richly-populated blue hori- 
zontal branch e s among Galactic globular clusters (e.g. , 
Menzied 1 19721: ISamus et al.l 1 19951 : [Melbourne et al.l 120001 : 
Piotto et al. 20021 1. This metal-poor clust er is a member of 
the “old” globular cluster population fe.g.. [Melbourne et al.l : 


* This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile, 
t Email: iur@umich.edu 
t Email: ian@obs.carnegiescience.edu 


iMarm-Franch et al.ll200^ 1 , and it ranks in the t op quartile o f 
most luminous clusters around the Milky Way llHarrislll996l l . 
The proper motion m e asurem ents and orbit calculations of 
ICasetti-Dinescu et al.l (l2007l l suggest that NGC 4833 and 
NGC 5986, another ol d and metal-poor cluste r, share similar 
orbital characteristics. ICasetti-Dinescu et al.l suggested that 
these two clusters are dynamically associated and could be 
an accreted pair with a common origin in a now-disrupted 
satellite galaxy. 


Until last year, the only high-resolution spectroscopic 
observations of stars in NGC 4833 h ad be en conducted by 
[Pilachowski, Sneden. fc Wallersteinl lll983li , who observed 


two stars, 
stars, and 


Gratton fc Ortolanil 


Minniti et al. 1199 


, who observed two 
, who observed one. 
15-20 species in indi- 


These studies derived abundances of 
vidual stars in NGC 4833, including O and Na abundances. 
They recognized that [O/Fe] in some stars was depleted rel¬ 
ative to the maximum [O/Fe] ratios found in other clusters 
and field stars of similar metallicity. The small sample sizes 
limited these authors’ abilities t o characterize the gl obal 
chemical properties of NGC 4833. [^rretta et al.l (l2014h ob¬ 
served 78 stars in NGC 4833 with the UVES and GIRAFFE 
spectrometers in the FLAMES instrument on the VLT. 
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Table 1. Basic properties and orbital parameters of NGC 4833 


Quantity 

Symbol 

Value 

Reference 

Right ascension 

a (J2000) 

12:59:34 

1 

Declination 

<5 (J2000) 

-70:52:35 

1 

Galactic longitude 

e 

303.6° 

1 

Galactic latitude 

b 

-8.0° 

1 

Luminosity 

Mv 

-8.17 

1 

Mass-to-light ratio 

M/L 

0.84 ± 0.45 

2 

Mass 

M* 

1.2x10® M© 

1 

Gentral concentration 

c 

1.25 

1 

Distance modulus 

m — M 

15.05 ± 0.06 

3 

Distance from Sun 

R© 

6.5 kpc 

4 

Distance from Galactic Center 

Rg.c. 

7.0 kpc 

4 

Perigalactic distance 

Rperi 

0.7 ± 0.2 kpc 

4 

Apogalactic distance 

Rapo 

7.7 ± 0.7 kpc 

4 

Maximum distance above Galactic plane 

Zmax 

1.8 ± 0.4 kpc 

4 

Orbital eccentricity 

e 

0.84 ± 0.03 

4 

Orbital period 

^orbit 

91 ± 9 Myr 

4 

References: (1) iHarrid Il996l: (2) ICarretta et alJ l2014l; 

(3) iMelbourne 

et al.l I 2 OOOI: 


This dataset enabled ICarretta et alJ to characterize the pat¬ 
tern of light-element (O, Na, Mg, Al, and Si) variations in 
NGC 4833. That study also characterized the abundances 
of 18 species of 15 heavier elements. 

ICarretta et al.l (l2014fl presented several forms of evi¬ 
dence that suggest that NGC 4833 may have been tidally 
stripped more than the average cluster. First, NGC 4833 
has a wide span of [Mg/Fe] ratios. This property has only 
been identif i ed in a few massive or me t al-poor clusters 
llYorig et al.l l2005l : ICarretta et al.l l2009bl : ICohen fc Kirbvl 
I2ni2h . vet NGC 4833 is only the 36*^ most luminous Galac¬ 
tic globular cluster today. Second, NGC 4833 has an ec¬ 
centric orbit, small perigalactic distance, and modest cen¬ 
tral concentration. These valu es predict a high destruction 
rate relative to other clusters (|Casetti-Dinescu et al.ll200n : 

lAllen. Moreno, fc Pichard9 120081 '). Third, the inner-guartile 

range of the [O/Na] ratios (IQR[0/Na]) in NGC 4833 lies on 
the upper envelope of clusters in the IQR[0 /Na1 versus lu¬ 
minos ity correlation. Evidence discussed by I Garrett a et al] 
(l2014h indicates that less-concentrated clusters of a given lu¬ 
minosity may have larger values of IQR[0/Na] and have lost 
more mass than more-concentrated ones. Finally, the low 
mass-to-light ratio of NGC 4833 hints that significant mass 
loss may have occurred as low- mass stars hav e been pref¬ 
erentially lost from the cluster. ICarretta et al.l (l2014ll note 
that no attempt has been made to detect these escaped stars 
as tidal tails, and potential investigators may have been in¬ 
timidated by the stellar crowding and high reddening at the 
low Galactic latitude of NGC 4833. The preponderance of 
evidence suggests that NGC 4833 has lost a larger fraction 
of its initial stellar mass than the average Galactic globular 
cluster. 

Our sample includes only 15 stars in NGC 4833, but 
our data cover most of the optical spectral range. We de¬ 
tect 43 species of 39 elements heavier than He, and this 
complements the more limited chemical in ventory derived 
from a larger sample of stars studied by ICarretta et all 
(l2014h . These data enable us to focus new attention on 
the abundance pattern exhibited by the heaviest elements 


in NGC 4833, those produced by neutron (n) capture reac¬ 
tions. 

Throughout this work we adopt the standard definitions 
of elemental abundances and ratios. For element X, the log¬ 
arithmic abundance is defined as the number of atoms of X 
per 10^^ hydrogen atoms, loge(X) = log^Q(Ax/AH)+ 12.0. 
For elements X and Y, [X/Y] is the logarithmic abundance 
ratio relative to the solar ratio, defined as logjQ(A x/Ay) — 
logj^g (Ax/Ay)©, using like ionization states; i.e., neutrals 
with neutrals a nd ions with ions. We adopt the solar abun¬ 
dances listed in lAsplund et al] (120091 '). Abundances or ratios 
denoted with the ionization state indicate the total elemen¬ 
tal abundance as derived from transitions of that particular 
state. 


2 OBSERVATIONS 

We have observed 15 probable members of NGC 4833 with 
the M agellan Inamori Kyoc era Echelle (MIKE) spectro¬ 
graph jBernstein et al.l l2003l i on the 6.5 m Landon Clay 
(Magellan II) Telescope at Las Campanas Observatory. Ta¬ 
ble [5] presents a log of these observations, and Eigure [T] 
highlights these 15 stars on a color-magnitude diagram of 
NGC 4833. ThAr comparison lamp spectra have been taken 
immediately preceding or following each observation. The 
red and blue arms of MIKE are split by a dichroic at 
« 4950A. This setup provides complete wavelength cover¬ 
age from 3350-9150A. Data reduction, extraction, and wave¬ 
length calibration have been performed using the CarPy 
MIKE data r eduction pipeline written by D. Kelson (see also 
lKelsorjr2003l L Continuum normalization and order stitching 
have been performed within the iraf environment. 

Six stars were observed with the (L.'TxSVO slit. This 
setup yields a resolving power of R ^ 41,400 in the blue and 
R ~ 36,300 in the red, as measured from isolated ThAr lines 
in the extracted comparison lamp spectra. Nine stars were 
observed during poorer seeing conditions with the l"0x5V0 
slit, which yields a resolving power of R ^ 30,500 in the blue 
and R ~ 25,900 in the red. Table [3] presents signal-to-noise 
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Table 3. Photometry, model atmosphere parameters, and S/N ratios 


Star_^ 

Star^ 

Star^ 

V 

(B - V)o 

O 

1 

Teff 

(K) 

log 9 
[cgs] 

Vt 

(km s“^) 

[M/H] 

S/N 

(3950 A) 

S/N 

(4550 A) 

S/N 

(5200 A) 

S/N 

(6750 A) 

2-185 

D102 


14.57 

0.85 

2.35 

4670 

1.55 

1.75 

-2.2 

25 

50 

50 

110 

2-277 

D109 


14.90 

0.87 

2.26 

4752 

1.72 

1.80 

-2.2 

25 

50 

50 

110 

2-882 

D193 


14.66 

0.84 

2.45 

4665 

1.57 

1.70 

-2.2 

20 

40 

50 

110 

2-918 

C49 

33027 

14.75 

0.83 

2.40 

4696 

1.63 

1.85 

-2.2 

20 

45 

40 

no 

2-1578 

C91 

31836 

14.12 

0.89 

2.26 

4595 

1.35 

1.75 

-2.2 

30 

50 

60 

140 

2-1664 

C92 


14.34 

0.86 

2.82 

4527 

1.32 

1.80 

-2.2 

35 

75 

75 

180 

3-742 

C129 

31370 

13.60 

1.05 

2.91 

4357 

0.93 

2.00 

-2.2 

25 

65 

75 

200 

3-772 

D167 


14.50 

0.89 

2.37 

4647 

1.50 

1.65 

-2.2 

10 

30 

40 

95 

3-1509 

C145 

31929 

14.38 

0.87 

2.52 

4591 

1.41 

1.85 

-2.2 

20 

45 

50 

120 

4-224 

D136 


13.51 

1.02 

2.70 

4385 

0.96 

2.05 

-2.2 

20 

45 

55 

120 

4-341 

D135 


13.56 

1.05 

2.69 

4397 

0.98 

1.95 

-2.2 

40 

95 

110 

250 

4-464 

D153 


14.48 

0.86 

2.45 

4628 

1.48 

1.90 

-2.2 

10 

40 

45 

115 

4-1255 

D49 


13.21 

1.14 

2.60 

4344 

0.84 

2.15 

-2.2 

35 

80 

90 

215 

4-1398 

C191 


14.23 

0.89 

2.56 

4556 

1.33 

1.85 

-2.2 

15 

45 

50 

115 

4-1706 

D41 


15.01 

0.78 

2.32 

4764 

1.76 

1.85 

-2.2 

25 

45 

50 

105 


^ This study, based on IPiotto et alJ ll2002l~) 

^ i MenziesI ill972l'l _ 

^ Carretta et al.l ll2014h 



Figure 1. The B — V versus V c o lour-m agnitude diagram us¬ 
ing photometry from IPiotto et alJ l|2002h . Stars observed with 
MIKE are marked with cir cles. The line is an isochrone from the 
PARSEC database (v. 1.2S: lBressan et~al .l2ni2l:IChen et al.l2niil 
computed for an age of 13 Gyr and metallicity Z = 0.00025. The 
isochrone is shifted to a distance modulus of 15.05 and a redden¬ 
ing of E{B — V) = 0.32. 


Table 2. Observing log 


Star 

Slit 

Date 


UTstart 

fexp 

(s) 

RV 

(km s~^) 

2-185 

T'O 

2011 

03 

18 

04:41 

2600 

211.0 

2-277 

T'O 

2011 

03 

18 

03:38 

3600 

201.7 

2-882 

I'.'O 

2011 

03 

18 

02:42 

3200 

209.8 

2-918 

I'.'O 

2011 

03 

18 

05:27 

3000 

209.8 

2-1578 

0'.'7 

2011 

03 

17 

02:52 

5400 

201.3 

2-1578 

0'.'7 

2011 

03 

17 

06:33 

1800 

201.2 

2-1664 

0"7 

2011 

03 

17 

04:30 

7200 

199.5 

3-742 

0'.'7 

2011 

01 

31 

07:43 

2400 

206.4 

3-772 

I'.'O 

2011 

03 

17 

23:40 

2800 

203.0 

3-1509 

I'.'O 

2011 

03 

18 

01:56 

2400 

207.4 

4-224 

0'.'7 

2011 

01 

28 

08:22 

1200 

201.5 

4-341 

0'.'7 

2011 

01 

28 

08:47 

2400 

202.2 

4-341 

0'.'7 

2011 

01 

29 

08:22 

2400 

202.5 

4-464 

I'.'O 

2011 

03 

18 

00:29 

2600 

297.6 

4-1255 

0'.'7 

2011 

01 

31 

08:29 

2400 

205.9 

4-1398 

I'.'O 

2011 

03 

18 

01:16 

2200 

200.1 

4-1706 

I'.'O 

2011 

03 

18 

06:20 

3400 

187.9 


(S/N) estimates based on Poisson statistics for the number 
of photons collected per pixel in the continuum. These values 
are measured from the co-added spectra. 

We adopt a star naming convention based on the 
IPiotto et al.l ll2002ll photometry files. This convention has 
the form “chip-ID.” The first value (“chip”) indicates the 
WFPC2 chip number, and the second value (“ID”) indicates 
the star identification number. 

Figure [2] illustrates a 25 A portion of the MIKE spectra 
of all 15 stars observed in our study. The spectra are sorted 
by increasing F-magnitude from top to bottom. Lines of 
nine species of six elements are marked. This illustrates the 
richness of blue spectra of red giant stars, even when the 
overall metallicity is less than 100 times the solar metallicity. 
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Figure 2. Small portions of the blue spectra of all 15 stars in our study. Identifications are given for 19 absorption lines. Velocity 
shifts have been removed, and the spectra have been shifted vertically for the sake of illustration. The stars are sorted by increasing 
V-magnitude from top to bottom. 


3 RADIAL VELOCITIES 

We measure the stellar radial velocity (RV) of each observa¬ 
tion with respect to the ThAr lamp by cross-correlating the 
echelle order containing the Mg l b lines in each spectrum 
against a template using the fxcor task in iraf. We create 
the template by measuring the wavelengths of unblended 
Fe I lines in this order in star 4-1255, which has the highest 
S/N in a single exposure. We compute velocity corrections 
to the Heliocentric rest frame using the iraf rvcorrect task. 
This method yie lds a total uncertaint y of « 0.6 km s“^ per 
observation (see iRoederer et al.ll20l4l . Table [5] presents the 
Heliocentric RV measured for each observation. 

We derive a mean Heliocentric RV to NGC 4833 of 
203.0 ± 1.5 km s“^ [s.e.m.] (cr = 5.9 km s“^ [s.d.]). All stars 
are within 2’ of the cluster center, and the RV of each star is 
within 2.5cr of the systemic velocity. There is a high proba¬ 
bility that all 15 stars observed are members of NGC 4833. 

The mean heliocentric RV and dispersion of NGC 4833 
have been derived from high resolution GIRAFFE spec- 
tra by two previo us studies , ICarre tta et all (l2014fl and 
iLardo et al.l (l2015t l. ICarretta et al.l reported a systemic 
velocity of 202.0 ± 0.5 km s”*^ and a d ispersion of 
4.1 ± 0.3 km s“^ based on a sample of 78 stars. Il^rdo et all 
reported a systemic velocity of 202.1 ± 0.6 km s“^ and a 
dispersion of 3.9 km s“^ based on a sample of 82 stars. 


Our values are in agreement. These values are also in 
reasonable agreement with previous measurements derived 
from low-resol ution spectra of sma ller samples of stars in 
NGC 4833 by iGeisler et al.l lll995h (201.0 rfc 1.3 km s~\ 
a — 4.6 kms“^, 12 stars) and iRutledee et al.l (Il997fl 
(194.1 ± 7.5 km s“^, a — 3.3 km s“^, 26 stars). For com¬ 
pleteness, we note that no ab undances hav e been presented 
from the spectra presented bv iLardo et al.l 


4 EQUIVALENT WIDTHS 

We measure equivalent widths using a semi-automatic rou¬ 
tine that fits Voigt absorption line profiles to the continuum- 
no rmalized spectra . Follo wing the same methods discussed 
in iRoederer et al.l ll2014l') . we inspect all equivalent width 
measurements by eye, and we make manual adjustments 
to the fits when necessary. Uncertainties in the equivalent 
widths are estimated based on the S/N ratios that vary as 
a function of wavelength. The complete list of equivalent 
widths is presented in Table 01 which is available in the Sup¬ 
plementary Information section found in the online edition 
of the journal. A short version is shown in the printed edition 
to illustrate its form and content. 

We o bserved four sta r s in co mmon with the GIRAFFE 
sample of ICarretta et al.l (l2014l l . Those authors note that 
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Table 4. Atomic data and equivalent widths 


Wavelength 

(A) 

Species^ 

E.P. 

(eV) 

logs/ 

Ref. 

2-185 ... 

6707.80 

3.0 

0.00 

0.17 

1 

limit 

6300.30 

8.0 

0.00 

-9.78 

2 

limit 

7771.94 

8.0 

9.14 

0.37 

2 

limit 

7774.17 

8.0 

9.14 

0.22 

2 

limit 

7775.39 

8.0 

9.14 

0.00 

2 

limit 

5682.63 

11.0 

2.10 

-0.71 

2 

23.9 

5688.20 

11.0 

2.10 

-0.45 

2 

43.3 

6154.22 

11.0 

2.10 

-1.55 

2 


6160.75 

11.0 

2.10 

-1.25 

2 

13.1 

4702.99 

12.0 

4.33 

-0.38 

2 

88.9 

5528.40 

12.0 

4.34 

-0.50 

2 

95.7 

5711.09 

12.0 

4.34 

-1.72 

2 

20.9 


Notes.— The full version of Tableis available in the supplemen¬ 
tary material online. Here, “synth” denotes lines used to derive 
an abundance via spectrum synthesis, and “limit” denotes lines 

used to derive an upper limit. _ 

Refer ences.— QJ _ [Smith. Lambert. NissenI 

19981: _ (2} _ iFuhr fc: Wiesel 120091: (3) 

Aldenius. Lundberg. fc Blackwell- Whitehead) _ [200^; (4) 


_ _ _ A) 

L^^e^^^Dakinl [li}89|, usin g HFS from Kunacz_&^en|_[l99^ 
(51 iLawleT'et'^ll l201.^ : (6}J|Woo^_et_alJ [20^; (7) 


2014 

a (9) 

Sobeck. Lawler. Kr. Snec 

enl I 2 OO 7 I: (lOi (Nilsson 

et al. 

2006 

; (11) 

Den Hartoe et al.l 2011 

for both log Q f value and 

HFS; (12) 

Ruffoni et al.l l2014l: ('lO^ 

iFuhr & Wiesel l2009l. us- 

ing 

HFS from iKurucz Belli Il995l 

(141 Wood et al.l 

2014al: 

iR.oederer Lawlerl l2012l: flB) 

Biemont et al.l 12011 

(17) 

Liune et alJ 

20061: ('18') iNilsson et al 

I 2 OIOI: (191 IPalmeri 

et al. 

2005 

; (20") iFuhr V Wiesel I 2 OO 9 I. using HFS/IS from iMcWilliam 

1998 

when 

available: (21) ILawler. 

Bonvallet. Sneden 

200l|. 

using HFS 

from Ilvans et al.l I 2 OO 6 I: (22) IRoederer et al. 

2 OI 1 I: 

{'23^ iLawler 

et al.l 2 OO 9 I; (2A) iLi et al.l l2007l. usine HFS 

from 

Sneden et alJ 120091: ('25') llvarsson. 

Litzen. WahlerenI 

2001 . 

usine HFS fromjSneden et alj 20091: f26)lDen Hartoe et al. 

2003. 


using HFS/IS fro m iRoederer et alJ [200^ when ava ilable; (27) 


iLawler et al.ll2006l. using HFS/IS from IRoederer et al. [2008j when 
avail able; (28 MLawler et al]l20Qla , usin ^ HFS/IS from Ivar^et al.l 
20061; (29') jPen Harto^^^al]l200t ; (30 HRoederer et al.|| 2012|: fSlI 


Lawler et al .1 l2n01bl. usings HFS from iLaw 'ler^^a^b^^ Blaisel 
200ll when available; (32) IWi ckliffe. Lawler. i^Na\^2000l; (33) 

Lawlerj_Snedenj_fc CowajJ [200^ for both log q f value and HFS; 

|200^ (35) IWickliffe fc Lawled Il997l : (36) 
Sneden^^^f^^M for b oth log q f val ue and HFS/IS; (37) 


Lawler et al.l l2007l : (381 llvarsson et alJ l2003l . using HFS/IS 
from IC owTiet alJ l2005l— see note on log of values there; (39) 
iBiemor^T^an' 20 00|, i7 ing HFS/IS from iRoederer et alJ 120121 : 
(401 InIIssou et al.ll2002l . 


“ The number to the left of the decimal point indicates the atomic 
number, and the number to the right of the decimal point indicates 
the ionization state (0 = neutral, 1 = first ion) 


equivalent widths measured from their GIRAFFE spectra 
are consistently larger than the equivalent widths measured 
from their higher-resolution UVES spectra. They adopt a 
linear transformation to convert their GIRAFFE equiva¬ 
lent widths to the UVES scale. We compare our measured 
equivalent widths with their corrected GIRAEEE equiva- 


T3 

3 


£ 

Cd 


»<; 

E 




Wavelength (A) 

Figure 3. Compa rison of equivalent widths measured by 
ICarretta et al.l ll2014h and us for four stars in common. The top 
panel illustrates the equivalent widths, and the bottom panel il¬ 
lustrates the differences as a function of wavelength. The dotted 
lines represent a 1:1 correspondence, and the solid lines represent 
the linear fits. 


lent widths. As illustrated in Figure O we find a difference of 
+4.7 ± 0.7 mA (a = 5.4 mA) in the sense that our equivalent 
widths are larger. There is a slight trend with wavelength, 
but this is only sign i ficant at the 2a level. 

ICarretta et al.l (l2007h note that the transformation 
from the GIRAFFE scale to the UVES scale is necessary to 
ensure agreement with equivalent widths measured indepen¬ 
dently by other investigators. Equivalent widths measured 
from their UVES spectra of stars in NGC 6752 are in rea¬ 
sonable agreement with those measu red from R ~ 110,000 
spectra obtained bv lYong et al.l (l2003l) . (A) = +1.7+0.4 mA 
(a = 5.6 mA). Our equivalent width measurements are 
made with the same techniques and software used to ex- 
amine hundreds of m etal-poor stars observed with MIKE 
(IRoederer et al.ll2014l) . The CarPy software package used 
to extract the MIKE observations is also identical. Our pre¬ 
vious work found excellent agreement among the equiva¬ 
lent widths measured from MIKE spectra and those mea¬ 
sured from echelle data collected at McDonald Observatory. 
There was a difference of +0.7 ± 0.2 mA {a — 4.0 mA) 
for the HRS on the Hobby-Eberly Telescope and a differ¬ 
ence of +1.0 ± 0.1 mA {a = 3.7 mA) for the Tull Coude 
Spect rograph on the Harlan J. Smith Telescope. Eurther- 
more, IRoederer et al.l made extensive comparisons between 
the equivalent widths f or weak lines in st ars in common wit h 
ICarretta et al.l ll2002h . Ijohnsoni (I 2 OO 2 I ). Ilvans etTP ll2003l) . 
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ICavrel et alj il2004l ~l. iHonda et alJ ll2004h . and iLai et all 
1 2 OO 8 II . The offsets were all ^ 1.3 mA. Thus we conclude 
that equivalent widths measured from MIKE spectra taken 
with the OV 7 slit are not systematically in error. 

We also investigate whether the wider MIKE slit used 
for some of our observations may be the source of this dif- 
ference. Two o f the fo ur stars in NGC 4833 in common with 
ICarretta et al.l (l2014h were observed using the 0”7 slit, and 
two were observed using the T.'O slit. The mean difference 
is +4.5 ± 1.1 mA {a = 6.2 mA) for the two stars observed 
with the 0”7 slit and +4.8 ± 0.8 mA (cr = 4.6 mA) for the 
two stars observed with the l"0 slit. The difference is not 
significant, and we conclude that the different slits are not 
the source of the discrepancy. 

In principle, poor sky subtraction could result in a sys¬ 
tematic bias in the equivalent widths measu red from a sin¬ 
gle ni ght. The four stars in common with I Garrett a et al.l 
ll2014ll were observed on three different nights. Star 3-742 
was observed on 2011 01 31, when the moon was only 6 per 
cent illuminated and 69° away from NGC 4833. The other 
three stars were observed on 2011 03 17 and 2011 03 18, 
when the moon was >90 per cent illuminated and >78° 
away from NGC 4833. The mean difference is +5.6±1.6 mA 
(cr = 6.7 mA) for the observations near new moon and 
+4.2 ± 0.7 mA (a = 4.6 mA) for the observations near full 
moon. The difference is not significant, and we conclude that 
poor sky subtraction is not the source of the discrepancy. 

_ Finally, we note that Roederer, Marino, & Sneden 

(l201iri examined spectra of several stars in globular clus¬ 
ter M22 that were taken on the same two nights in March, 
2011. That study compared the equivalent widths mea¬ 
sured in those stars with those measured by iMarino et al.l 
(l201llf . Again, the offset is small, (A) = +1.8 ± 0.3 mA 
(a = 3.7 mA). There is no compelling evidence that our 
MIKE observations of M22 or NGC 4833 conducted in 
March, 2011, are contaminated by the solar spectrum re¬ 
flected by the full moon. 

In conclusion, we cannot explain the 4.7 mA offset in 
the me an equivalent width di fferences between our data and 
that of I Garrett a et all (l2014lf . In Section [6.21 we find that 
the cumulative effect of the choice of Fe i lines used for 
the analysis, the log gf values, and this equivalent width 
difference is 0.04 dex. This value is relatively small, and we 
do not pursue the matter further. 


5 BROADBAND PHOTOMETRY AND 
REDDENING 

TableElpresent s the photome t ry an d colours. V and (B — V) 
are taken from IPiotto et al.l (l2002l 'n and have been trans¬ 
formed by those authors to the Johnson system from their 
Hubble Space Telescope (HST) WFPC2 FA39W and F555W 
broadband photometry. K magnitu des are taken from th e 
Two-Micron All Sky Survey l2MASS: ISkrutskie et al.ll2006ll . 
The {B — V) and {V — K) colours a re dereddened according 
to the differential reddening maps of[Melboumeet^iJj 20001 1 
and the extinction coefficients of ICardelli et akF i 1981 ). All 


^ Downloaded from http://www.astro.unipd.it/globulars/ on 
2010 October 22 



(V-K)„ 

3.11 

2.82 

2.53 

2.24 

1.95 


Figure 4. Relationship between dereddened K magnitudes and 
Teff predicted by dereddened {V — K) colours, which are approx¬ 
imated on the right axis for reference. The linear least-squares fit 
is shown. 


15 stars that we have observed are within 2’ of the clus¬ 
ter center, and the full range of reddening for these stars is 
0.310 ^ E{B — V) ^ 0.322 according to the |Melbourne et al] 
maps. The differential reddening is a small effect compared 
to the relatively large extinction toward NGC 4833 at low 
Galactic latitude (fe = —8°). 


6 MODEL ATMOSPHERES 

We interpolate model atmospheres fro m the ATLAS9 a- 
enhanced grid llGastelli fc Kuruc3l2003l) . We use a re cent 
versi on of the line anal ysis code MOOG llSnedenlll973l : see 
also ISobeck et aDl201lf) . operated in batch mode, to per¬ 
form the analysis. In Section [6.II we discuss our methods to 
estimate the appropriate model atmosphere parameters. We 
compa re these values with those derived bv iGarretta et ahl 
ll2014ll in Section 


6.1 Model parameters 

We first estimate the stellar effective temperature (Teft) us- 
ing the metal l icity-d ependent V — K colour-Tefi relations of 
lAlonso et al.l (Il999a) . We transform the K magnitudes from 
the 2MA SS system to the TCS syste m according to Equa¬ 
tion 5 of iRamfrez fc Melended ll2005ll . We refine the initial 
estimate by fitting a line to the relation between dereddened 
K magnitude, Kq, and Teff calculated from {V — K)q. Fig¬ 
ure |4] illustrates this relationship. The standard deviation of 
the residuals with respect to this fit is 97 K, which we adopt 
as the uncertainty in TeS for each star. 

The two stars that deviate most from the fit are 2-1578 
and 2-1664. Both stars are within 2' of the cluster cen¬ 
ter, as are most other stars in our sample. No neighboring 
stars within 5 mag nitudes and 10 pixels are found in the 
IPiotto et al.1 ll2002l ') photometry files, and no faint contami¬ 
nants are apparent in the WFPC2 images. No photometry 
error flags are reported for either star. The {V — K)o colors 
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Abundances in Globular Cluster NGC 4833 7 


of these stars would need to be different by about +0.2 mag 
and —0.2 mag, respectively, to bring each star to the best-fit 
line in Figure [3] These stars are only about 11” apart, so 
it is highly unlikely that differential reddening is underesti¬ 
mated for one and overestimated for the other. We derive 
metallicities of these two stars ([Fe l/H] = —2.23 and —2.24, 
[Fe ll/H] = —2.23 and —2.18, respectively; Section that 
are in excellent agreement with the cluster means. In sum¬ 
mary, we find no reason to exclude either star when comput¬ 
ing the standard deviation of the residuals to the fit shown 
in Figure |4] 

Uncertainties in Kq translate to uncertainties of only 
33 K in Teff on average. Small changes in the differen¬ 
tial reddening across the face of NGC 4833 have a min¬ 
imal effect on when derive d this way. A{K)/ E{B — 
V) « 0.35 (ICardelli et al.|[l989ll . so a differential redden¬ 
ing of ± 0.012 mag in E{B — V) dMelbourne et al.l bood l 
translates to ± 0.004 mag in K. The slope of the Tes-Ko 
relation shown in Figure [4] is 201 K/mag. Thus differential 
reddening contributes an average of « 1 K to the Tefr uncer¬ 
tainty budget, which is negligible. We assume that all giants 
in NGC 4833 reside along a single red giant branch (RGB) 
fiducial with [Fe/H] = —2.2 ('Section lO.lll . A change in [Fe/H] 
of 0.1 dex affects the calculated Tefj values by < 1 K, which 
is al so negligible. _ 

iMelbourne et all (120001 1 derived the distance to 
NGC 4833 using the absolute magnitudes of RR Lyrae vari¬ 
ables along the horizontal branch (HB). We use this dis¬ 
tance, (m — M) = 15.05 ± 0.06, to derive the surface grav¬ 
ities (log g) for all 15 stars observed from the fundamental 
relation logg = 41og(Teff,*) - log(M/M©) - q - 

Mhoi,*) — 10.61. Here, Mboi,* = BCic + niK — {m — M), and 
the consta nt 10.61 is derived from the solar values given in 
ICoxI (!2000I 1. We assume a mass of 0.8 M© for all stars on the 
RGB in NGC 4833. We interpolate BCrc from the grid o f 
bolometric corrections presented bv lAlonso et al.l (Il999al '). 
Uncertainties of 0.05 in magnitude, 97 K in Tefr, 0.1 M© 
in mass, and 0.06 in the distance modulus only affect the 
derived log g values by 0.02, 0.04, 0.05, and 0.03 dex, re¬ 
spectively. Their collective uncertainties contribute a total 
of 0.07 dex uncertainty in log g. 

We derive the microturbulent velocity, ut, by minimiz¬ 
ing correlations between the equivalent width (expressed as 
the log of the equivalent width divided by A) of Fe l lines and 
the abundances derived from them. The internal uncertainty 
associated with this method is « 0.05 dex. 

Finally, we assign a single metallicity to all model at¬ 
mospheres, [Fe/H] = —2.2. This is justified by our results 
indicating that NGC 4833 is mono-metallic and has a mean 
metallicity of [Fe ii/H] = —2.19 ±0.013 fSection 19.11) . Our 
use of the a-enhanced grid of models is justified because the 
most abundant electron-donating metals are enhanced by 
factors of a few relative to Fe in NGC 4833 when compared 
with the solar ratios. 


6.2 Comparison with Carretta et al. 

ICarretta et al.l (l2014h used a procedure similar to ours to 
derive their Tefr values for sta rs in NGC 4833 . Thei r proce¬ 
dure is discussed in detail by iGratton et al.l ll2007lf . There 
are four stars in common between our study and theirs, and 
we derive a mean difference in Tefi of —28 ± 1 K (u = 2 K). 


Both studies use the same 2MASS K magnitudes to cal¬ 
culate Teff after computing a relation like that shown in 
Figure 01 There appears to be a zeropoint offset in the V 
magnitudes used to compute the initial relation ship between 
V — K and Teff. We adopt V magnitudes from IPiotto et aP 
1 I 2 OO 2 II . who calibrated their HST F555W magni tudes to V 
in the standard Johnson system. ICarretta et al.l used their 
own ground-based V magnitudes calibrated with an absolute 
zeropoint uncertainty of 0.03 mag. A systematic zeropoint 
difference of ± 0.03 mag translates to an average difference 
in Teff of + 27 K, exactly the difference foun d for these four 
stars. The Tgff offset between our study and ICarretta et ^ 
can be explained by the stated uncertainties in the U-band 
zeropoints. 

On avera ge, our log q values are lower than those 
calculated by TCarretta et al.l (l2014lf by 0.21 ± 0.02 dex 
(cr = 0.04 dex). Small differences in the adopted solar con¬ 
stants, mass of stars on the RGB in NGG 4833, and distance 
modulus to N GG 4833 only affect this offset by 0.04 dex. 
Adopting the iGarretta et al.l Teff values changes the log g 
values by only 0.01 dex. We adopt 0.21 dex as the system¬ 
atic uncertainty in log g. 

Our vt values are greater by 0.18 ± 0.04 km s“^ 
(fj = 0.08 km s~^) for the four stars in common with 
IGarretta et al.l (l2014l l. We adopt 0.18 km s“^ as the sys¬ 
tematic uncertainty in vt- 

Our metallicities (derived from Fe l lines) are lower by 
0.20 ± 0. 02 dex (a = 0. 0 4 dex ) for the four stars in com¬ 
mon with ICarretta et al.l (120141 '). The stated offset accounts 
for the different adopted solar Fe abundances. The discrep¬ 
ancy can be attributed to the model atmosphere grids and 
line analysis codes, as we now sh ow. We rederive metallici¬ 
ties for these four stars using the IGarretta et al.l model at¬ 
mosphere parameters, their line list and equivalent widths, 
our adopted grid of models (atlas9), and our adopted 
line analysis code (moog). The metallicities are lower by 
0.20 ± 0.01 dex (ct = 0.02 dex), which fully accounts for the 
discrepancy. 


7 ABUNDANCE ANALYSIS 


Abundances of lighter elements {Z ^ 30) that do not ex¬ 
hibit broad isotope shifts (IS) or hyperfine splitting structure 
(HFS) are derived from equivalent widths. We use MOOG 
to compute theoretical equivalent widths that are forced to 
match the measured values by adjusting the abundance. 
Abundances of all other elements are derived by match¬ 
ing synthetic spectra with the observed spectra. Lines of 
these elements are denoted by the word “synth” in Table 01 
When a line is not detected, we derive a 3a abundance up- 
pe r limit using a versio n of the formula presented on p. 590 
I), which is derived from equation A8 
These lines are denoted by the word 


Frebel et al.l 


Bohlin et al.l 


H 


of 
of 

“limit” in Table 

The atomic data for all lines examined are presented 
in Table 01 References for the loggf values, IS, and HFS 
(when available) are also given in Table 01 We derive C 
abundances from the CH A^A — X^Il G band (« 4290- 
4330 A) using a line list kindly provided by B. Plez (2007, 
private communication). We derive N abundances from the 
GN — X^E violet band (« 3875-3885 A) using the line 
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list from iKurucz lll995ll after setting the C abun¬ 

dance using the CH G band. 

We assume local thermodynamic equilibrium (LTE) 
holds for the line-forming layers of the atmosphere for all 
species except for tho se described below . We adopt non-LTE 
corre ctions for Li l jLind e t al.l 1200^ . 0_j_j[Fabbjanet^lJ 
[ 200 ^, Na I (iLind et al.ll201ll l. and K l llTakeda et al.l 2002ll . 
We use the correction for the nearest point on the grid when 
the stellar parameters extend beyond the edges of the grids 
of pre-computed non-LTE corrections. These corrections are 
included in the abundances presented in Tables O |6] and [3 
We include an additional 0.1 dex statistical uncertainty in 
our error estimates to account for uncertainties in the cor¬ 
rections. We do not detect the Li l 6707 A line; corrections 
are based on the Scr upper limits and range from -fO.lO to 
-1-0.21 dex. The non-LTE corrections for each of the three 
high-excitation O l 7771, 7774, and 7775 A triplet lines range 
from —0.05 to —0.08 dex. The non-LTE corrections for the 
Na I 5682, 5688, 6154, and 6160 A lines range from —0.02 
to —0.14 dex. The K l 7698 A line is contaminated by an 
atmospheric O 2 line; fort unately, however, the K l 7664 A 
line is not contaminated. iRoederer et al.l (|20l3) found that 
these two lines yield consistent abundances (to better than 
0.02 dex) after non-LTE corrections are applied. The non- 
LTE corrections for the K l 7664 A line range from —0.58 
to —0.61 dex. 

We apply one additional correction to the O abun¬ 
dance derived from the high-excitation O l triplet lines. The 
[O l] 6300 A line is generally considered to be a reli able abun¬ 
danc e indicator formed under conditions of LTE (iKiselmanI 
Hool), but the O i 7771, 7774, and 7775 A lines are not. 
Both abundance indicators are detected in seven stars. Af¬ 
ter applying the non-LTE corrections for the O l triplet lines, 
these lines still yield abundances higher by -1-0.40 ± 0.05 dex 
((7 = 0.13 dex) than the [O l] 6300 A line. This offset be- 
tw een the O l and [O l ] lines is c omparable to those fou nd 
bv iGarcfa Perez et al.l ( 20061 ) and IRoederer et ^ ll2014l ~l in 
other metal-poor giants. We apply a downward correction 
to all non-LTE-corrected O l triplet line abundances by 
0.40 dex. We include an additional 0.05 dex uncertainty to 
account for the uncertainty in this correction. 


8 RESULTS 

Table [5] lists the abundances or upper limits derived from 
each line examined in each star. Table [6] lists the weighted 
mean abundances of each species examined in each star. 
We adopt the upper limit that provides the strongest con¬ 
straint on the abundance when multiple lines of the same 
species are not detected. The complete versions of Tables [S] 
and [6] are available in the Supplementary Information avail¬ 
able online. Table [7] lists the mean abundance ratios of all 
species examined in NGG 4 833. These are computed using 
the formalism presented in iMcWilliam et all (Il995ll . The 
statistical uncerta inty, (Tstatisticai, is given by equation A17 of 
iMcWilliam et aP This includes uncertainties in the equiva¬ 
lent width measurement or line profile fitting, log gf values, 
and any non-LTE correctio ns. The total unce rtainty, ctotai, 
is given by equation A16 of lMcWilliam et al.l This includes 
the statistical uncertainty and uncertainties in the model 
atmosphere parameters. We recommend that (Jneutrais for 


Table 5. Abundances derived from individual lines 


Star 

Species 

Wavelength 

(A) 

log(E 

(T 

2-185 

Li I 

6707.80 

< 0.53 


2-185 

[O I] 

6300.30 

< 7.16 


2-185 

0 I 

7771.94 

< 6.83 


2-185 

0 I 

7774.17 

< 7.24 


2-185 

0 I 

7775.39 

< 7.40 


2-185 

Na I 

5682.63 

4.51 

0.20 

2-185 

Na I 

5688.20 

4.60 

0.20 

2-185 

Na I 

6160.75 

4.74 

0.21 

2-185 

Mg I 

4702.99 

5.38 

0.32 

2-185 

Mg I 

5528.40 

5.60 

0.29 

2-185 

Mg I 

5711.09 

5.53 

0.15 


Notes.— The full version of Table 0 is available in 
the Supplementary Information section online. 


element A be added in quadrature with o-gtatisticai for ele¬ 
ment B when computing the ratio [A/B] when B is derived 
from lines of the neutral species. We recommend a similar 
procedure, utilizing Uions instead of (Tneutrais, when element 
B is derived from lines of the ionized species. The adopted 
solar reference abundances are listed in Table El 

Figures[5]and[6]illustrate abundance ratios as a function 
of Teff. In most cases, no trends are found. Subtle trends may 
be present in cases where one of the elements in question 
is known to vary with evolutionary status. A small down¬ 
ward trend of [G/Fe] and a corresponding increase in [N/Fe] 
are found with decreasing Teff. The Li upper limits get pro¬ 
gressively lower with decreasing Tea, which is simply a re¬ 
flection of relatively constant equivalent width upper limits 
(from relatively constant S/N levels) applied to stars with 
decreasing Tea- There are also subtle differences at the « 1- 
2(7 level (« 0.10-0.15 dex) in the mean [X/Fe] ratios between 
the warm {Tea > 4500 K) and cool {Tea < 4500 K) stars, 
where “X” represents Sc, V l and ll, Zn, Nd, Ce, Sm, Eu, and 
Dy. We investigate these differences in detail in Sections 19.31 
and 19.51 

Table E lists the variations in abundance ratios that 
result from changes in the model atmosphere parameters 
for one representative star with high S/N ratios, 4-1255. The 
final column of Table [9] lists the total variation, calculated 
by adding the four individual variations in quadrature. 

8.1 Comparison of abundance ratios 

Table [3 also l i sts th e mean abundance ratios derived by 
iGarretta et al.l (l2014h from their UVES spectra of stars in 
NGG 4833. In many cases the two sets of mean abundances 
are in reasonable (2(7) agreement. At first glance, several 
[X/Fe] ratios are not in agreement, where “X” represents 
O I, Si I, Ca I, Sc II, Ti l and ll, V l, Gr ll, Co l, Ni l, Zn l, 
and Nd ll. Most of these differences can be explained by the 
sets of lines used, the loggf values, and the adopted solar 
abu ndances, as we n ow sh ow. 

ICarretta et al.l |20l3) used the atomic data and so - 
lar reference abundances presented in iGratton et al.l (l2003ll . 
Their approach ensures a uniform abundance scale for the 
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Table 6. Mean abundances in each star 


Star 

Species 

-^lines 

loge 

[X/Fe]^ 

^statistical 

•^total 

'^neutrals 

*^ions 

2-185 

Fe I 

102 

5.33 

-2.17 

0.06 

0.14 

0.00 

0.00 

2-185 

Fe II 

6 

5.27 

-2.23 

0.07 

0.10 

0.00 

0.00 

2-185 

Li I 

1 

< 0.53 






2-185 

C (CH) 

1 

5.62 

-0.64 

0.10 

0.22 

0.15 

0.15 

2-185 

N (CN) 

1 

7.29 

1.63 

0.20 

0.28 

0.23 

0.23 

2-185 

O I 

4 

< 6.83 

0.31 





2-185 

Na I 

3 

4.61 

0.55 

0.12 

0.17 

0.13 

0.17 

2-185 

Mg I 

3 

5.52 

0.09 

0.07 

0.14 

0.09 

0.15 


Notes.— The full version of Table|6]is available in the Supplementary Information section 
online. 


“ [Fe/H] is given for Fe I and Fe II 



4800 4600 4400 4200 


T.„ (K) 




4800 4600 4400 4200 

T.„ (K) 



Figure 5. Derived [X/Fe] ratios (loge for Li) as a function of Tefi for Li through Cr I. The dotted lines represent the solar ratios, and 
the downward facing triangles represent upper limits. 


dozens of globular clusters observed in the last decade using 
GIRAFFE and UVES. Our analysis incorporates atomic 
data from mo r e rece nt laboratory studies and adopts the 
lAsplund et al.l (l2009l) solar abundance scale. The appeal of 
using updated log gf values is that recent laboratory studies 
regularly report uncertainties on absolute transition proba¬ 
bilities as small as 5 per cent, which is effectively negligi¬ 
ble when compared with other sources of uncertainty. The 
choice of loggr/ values is irrelevant if only relative differ¬ 
ences are of interest and if the same set of lines is employed 
for each star. This condition is not met since the S/N and 
stellar parameters vary, thus we use updated loggf values 
whenever possible. 


The [Ca/Fe] ratios provide a good example of the 
consequences of this choice. There are seven Ca l 
lines studied i n the GIRAFFE spectra obtained by 
iGarretta et al.l (l2014l ). The Ca i log gf values used 

by CarTetta^_eL_alJ and those we adopt fr om _ NIST or 

lAldenius. Lundberg. fc Blackwell- Whiteheadl ll2009h differ 


by —0.26 dex to -1-0.27 dex for these lines. Only one of these 
absolute differences is smaller than 0.17 dex. In our study, 
we use three of these lines and several others at shorter 
wavelengths. Thus the exact set of lines measured signif¬ 
icantly affects the derived [Ca/Fe] ratios. Substantial dif¬ 
ferences 0.05 dex) are also present in the adopted log gf 
values of lines of Na i, Si i, Sc ii, Ti i, Fe ii, and Ni i. It is im¬ 
possible to assess the exact differences since no single source 
of log gf values covers all of the lines used bv ICarretta et aP 
and us. 

Our use of the lAsplund et al.l (l2009l ) solar abundance 
scale introduces s everal substanti a l diffe rences in the [X/Fe] 
ratios relative to ICarretta et al.l ll20l4l . Our [O i/Fe] ra¬ 
tios would change by —0.06 dex on their scale; [Na i/Fe], 
+0.07 dex; [Mg l/Fe], +0.21 dex; [Ca l/Fe], +0.11 dex; 
[Ti n/Fe[, —0.13 dex; [Cr ll/Fe], —0.08 dex; [Mn l/Fe], 
+0.13 dex; and [Zn i/Fe], +0.07 dex. 

We use our software, model atmospheres, and solar ref¬ 
erence abundances to rederive the [X/Fe] ratios using the 
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Table 7. Mean abundance ratios in NGC 4833 


Ratio 

Species or 
molecule 


This study 



Carretta et al. {'20141 

Mean 

Std. err. 

Std. dev. 

Num. 

stars 

Notes 

Mean 

Std. err. 

[C/Fe] 

GH 

-0.44 

0.06 

0.24 

15 




[N/Fe] 

GN 

+1.28 

0.09 

0.36 

15 




[O/Fe] 


+0.81 

0.08 

0.22 

8 

1 

+0.17 

0.22 

[Na/Fe] 


+0.39 

0.06 

0.22 

14 

1 

+0.52 

0.29 

[Mg/Fe] 


+0.35 

0.06 

0.24 

15 


+0.27 

0.22 

[Al/Fe] 


+1.09 

0.07 

0.20 

8 

1 

+0.90 

0.34 

[Si/Fe] 


+0.74 

0.02 

0.09 

15 


+0.47 

0.04 

[K/Fe] 


+0.40 

0.02 

0.10 

15 




[Ca/Fe] 


+0.49 

0.011 

0.04 

15 


+0.35 

0.01 

[Sc/Fe] 

II 

+0.19 

0.03 

0.11 

15 


-0.11 

0.01 

[Ti/Fe] 


+0.08 

0.014 

0.06 

15 


+0.18 

0.02 

[Ti/Fe] 

II 

+0.39 

0.02 

0.08 

15 


+0.23 

0.01 

[V/Fe] 


-0.44 

0.03 

0.10 

15 


-0.08 

0.01 

[V/Fe] 

II 

+0.05 

0.05 

0.19 

14 




[Cr/Fe] 


-0.26 

0.02 

0.07 

15 


-0.24 

0.02 

[Cr/Fe] 

II 

+0.29 

0.03 

0.12 

12 


+0.01 

0.05 

[Mn/Fe] 


-0.58 

0.03 

0.11 

15 


-0.54 

0.01 

[Fe/H] 


-2.25 

0.02 

0.09 

15 


-2.02 

0.01 

[Fe/H] 

II 

-2.19 

0.013 

0.05 

15 


-2.01 

0.02 

[Co/Fe] 


-0.39 

0.04 

0.14 

14 

1 

-0.03 

0.03 

[Ni/Fe] 


-0.06 

0.015 

0.06 

15 


-0.18 

0.01 

[Cu/Fe] 


-0.65 

0.05 

0.14 

9 

1 

-0.80 

0.09 

[Zn/Fe] 


+0.19 

0.02 

0.10 

15 


+0.07 

0.03 

[Rb/Fe] 


< +0.85 



15 




[Sr/Fe] 

II 

+0.18 

0.03 

0.12 

14 




[Y/Fe] 

II 

-0.22 

0.03 

0.11 

15 


-0.15 

0.06 

[Zr/Fe] 

II 

+0.25 

0.03 

0.11 

14 




[Nb/Fe] 

II 

< +0.93 



2 




[Mo/Fe] 


+0.19 

0.11 

0.23 

4 

1 



[Ba/Fe] 

II 

-0.02 

0.016 

0.06 

15 


-0.06 

0.07 

[La/Fe] 

II 

+0.08 

0.03 

0.13 

14 


+0.05 

0.03 

[Ce/Fe] 

II 

+0.03 

0.03 

0.11 

15 




[Pr/Fe] 

II 

+0.19 

0.09 

0.29 

10 

1 



[Nd/Fe] 

II 

+0.18 

0.03 

0.11 

15 




[Sm/Fe] 

II 

+0.29 

0.03 

0.11 

14 




[Eu/Fe] 

II 

+0.36 

0.03 

0.13 

15 




[Gd/Fe] 

II 

+0.32 

0.09 

0.21 

6 

1 



[Tb/Fe] 

II 

+0.19 

0.04 

0.05 

2 

1 



[Dy/Fe] 

II 

+0.43 

0.05 

0.17 

13 

1 



[Ho/Fe] 

II 

+0.16 

0.18 

0.18 

1 

1 



[Er/Fe] 

II 

+0.33 

0.08 

0.20 

7 

1 



[Tm/Fe] 

II 

+0.19 

0.07 

0.10 

2 

1 



[Yb/Fe] 

II 

+0.01 

0.08 

0.25 

9 

1 



[Hf/Fe] 

II 

+0.27 

0.04 

0.06 

2 

1 



[Ir/Fe] 


+0.66 

0.11 

0.11 

1 

1 



[Pb/Fe] 


< +0.02 



15 




[Th/Fe] 

II 

< +0.32 



14 





Notes— 1. Does not include one or more upper limits 


ICarretta et alj ll20l4) line lists and equ ivalent widths fo r the 
four stars in common. On our scale, the lCarretta et ahl ratios 
change by +0.15 ± 0.04 dex for [O i/Fe], +0.10 ± 0.03 dex 
for [Si i/Fe], —0.05 ± 0.05 dex for [Ca i/Fe], +0.05 ± 0.03 dex 
for [Ti l/Fe], and +0.07 ± 0.03 dex for [Ni l/Fe]. These ac¬ 
count for about half of the total discrepancy in O l, Si l, 
and Nil, but they widen the discrepancy for Ca l and Ti l. 
V I a nd Co l cannot be assessed due to the corrections for 
HFS. ICarretta et al.l onlv derived abundances for Ti ll, Cr ll. 


Zn I, and Nd ll from their UVES spectra, so we have no data 
in common. 

To summarize, some mean [X/Fe] rati os we have de¬ 
rived are in disagreement with those found bv ICarretta et all 
ll2014ll . We can explain some of these discrepancies in part 
by the different solar reference abundances, and we attribute 
others to the different sets of lines used and the log gf values 
of those lines. Finally, the exact sets of stars observed could 
affect the mean values of the light elements known to vary 
star-by-star. 
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Figure 6 . Derived [X/Fe] ratios ([Fe/H] for Fe) as a function of for Cr II through Yb. The dotted lines represent the solar ratios, 
and the downward facing triangles represent upper limits. 


9 DISCUSSION 
9.1 Iron 

We derive a mean metallicity, ([Fe/H]), of —2.25 ± 0.02 
(cr = 0.09) from Fe l lines in the 15 stars observed in 
NGC 4833. We derive ([Fe/H]) = -2.19 ± 0.013 (a = 0.05) 
from Fe ll lines. The standard deviation associated with each 
of these values does not exceed that expected from the model 
atmosphere uncertainties listed in Table |9] We Hnd no evi¬ 
dence for an intrinsic dispersion in the Fe abundances from 
one star to another within NG C 4833, which reaffirms the 
results of ICarretta et all ll2014h . 

As noted in Section IQ we find a meta llicity lower by 
0.20 ± 0.02 dex (cr = 0.04 dex) relative to ICarretta et al.l 
ll2014ll for four stars in common. We have identified the 
causes of this discrepancy previously, and now we compare 
the metallicity of NGC 4833 to an external metallicity stan¬ 
dard. The temperature of the K-giant Arcturus (a Boo) is 
known to better than 30 K from measurements of its angu - 
lar diameter and bolometric flux. lKoch fc McWilliamI (l2008l l 
performed a differential abundance analysis between Arc¬ 
turus and red giants in the moderately metal-poor globular 
cluster 47 Tu c ([Fe/H] = —0.76 ± 0.01 [statistical] ± 0.04 
[systematic! 1. liCrch fc McWilliarnl (1201 ill extended this dif¬ 
ferential analysis to the very metal-poor cluster NGC 6397 
([Fe/H] = —2.10 ± 0.02 [statistical] ± 0.07 [systematic]), 
which we use for our comparison. 

W e use the Fe l line list from iKoch fc McWilliamI 
1 I 2 OIII) to identify 19 lines in common to three red giants 
in NGC 6397 (stars 7230, 8958, and 13414) and three red 
giants in NGC 4833 (stars 3-742, 4-341, 4-1255) with sim¬ 
ilar Tefl. We measure equivalent widths for each of these 
lines in NGC 4833, and we perform a differential abun¬ 
dance analysis between NGC 6397 and NGC 4833. We hnd 
a mean difference in [Fe/H] of -1-0.09 ± 0.04 dex be¬ 
tween NGC 6397 and NGC 4833, In other words, if we 
adopt [Fe/H] = —2.10 ± 0.02 for NGC 6397 as derived by 


iKoch fc McWilliamI (1201 il l, we infer [Fe/H] = —2.19 ± 0.04 
for NGC 4833. This value is in fair agreement with our de¬ 
rived mean metallicity, [Fe/H] = —2.25 ± 0.02. On the basis 
of the di fferential globular c luster abundance scale estab¬ 
lished bv iKoch fc McWilliamI relative to Arcturus, we prefer 
the lower metallicity value for NGC 4833. 


9.2 Lithium through Silicon 

We do not detect Li in any of the stars in our sample. The 
upper limits we derive are illustrated in Figure [5] These 
limits rule out substantial enhancement of the surface Li 
abundances in these stars, as has been f ound in a few red 
giants in other globular cluste rs (e.g., iKraft et al.l Il999l : 
ISmith. Shetrone. fc Kean3ll999l ). 

Figure [7] shows the relationships among the light ele¬ 
ments C, N, and Na in NGC 4833. There is a clear decrease 
in the [C/Fe] ratio with increasing [N/Fe]. [C/Fe] also anti¬ 
correlates with [Na/Fe], and [N/Fe] correlates with [Na/Fe]. 
There is no compelling evidence for a bi-modality in either 
[C/Fe] or [N/Fe], but this tentative conclusion should be 
checked using larger sample sizes. 

Figure [S] reveals a prominent anti-correlation between 
[0/Fe] and [Na/Fe] in our sample of 15 stars in NGC 4833. 
ICarretta et al.l (l2014ll detected O and Na in 51 stars in 
their sample, and these data are shown for comparison in 
Figure El The offset between our mean [O/Fe] and that of 
ICarretta et al.l is apparent. This offset is strongly influenced 
by two factors. First, we only deriv e upper limits fo r [O/Fe] 
in stars with [O/Fe] < -1-0.4, while ICarretta et al.l detected 
O in « 20 such stars. These upper limits are not included 
in the mean [O/Fe] value reported in Table [T] Second, two 
of the eight stars with O detections in our sample show 
[O/Fe] « -1-1.0, which is « 0.2 dex higher than the next- 
highest [O/Fe] ratio in any other star in our sample. Nei- 
ther of these s tars, 4-224 and 4-1225, was included in the 
ICarretta et al.l sample. The remaining stars in our sample 
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12 Ian U. Roederer and Ian B. Thompson 


Table 9. Variations in Abundance Ratios Resulting from Changes in the Model 
Atmosphere Parameters for Star 4-1255 


Ratio 

Ion 

S[X/Y]/5nff 

<5[N/Y]/51og g 

S[X/Y]/5vt 

<5[X/Y]/(S[M/H] 


[O/Fe] 

I 

-0.27 

+0.10 

-1-0.06 

-1-0.04 

0.30 

[Na/Fe] 

I 

-0.06 

+0.00 

-1-0.07 

-1-0.01 

0.09 

[Mg/Fe] 

I 

-0.06 

-0.03 

-1-0.01 

-1-0.00 

0.07 

[Al/Fe] 

I 

-0.08 

-1-0.01 

-1-0.07 

-1-0.01 

0.11 

[Si/Fe] 

I 

-0.11 

-1-0.04 

-1-0.07 

-1-0.03 

0.14 

[K/Fe] 

I 

-1-0.03 

-1-0.01 

-0.03 

-0.01 

0.04 

[Ca/Fe] 

I 

-0.01 

-1-0.02 

-0.03 

-0.01 

0.04 

[Sc/Fe] 

II 

-1-0.01 

-1-0.05 

-1-0.07 

-1-0.05 

0.10 

[Ti/Fe] 

I 

-1-0.08 

-0.01 

-1-0.02 

-0.02 

0.09 

[Ti/Fe] 

II 

-0.06 

-0.02 

-0.04 

-0.01 

0.08 

[V/Fe] 

I 

-1-0.11 

-1-0.02 

-1-0.07 

-1-0.05 

0.14 

[V/Fe] 

II 

-0.02 

-0.05 

-1-0.00 

-0.07 

0.09 

[Cr/Fe] 

I 

-1-0.04 

-0.01 

-1-0.02 

-0.01 

0.05 

[Cr/Fe] 

II 

-0.12 

-1-0.00 

-1-0.03 

-0.02 

0.13 

[Mn/Fe] 

I 

-1-0.03 

-1-0.04 

-1-0.10 

-1-0.05 

0.12 

[Fe/H] 

I 

-1-0.15 

-0.03 

-0.07 

-0.02 

0.17 

[Fe/H] 

II 

-0.06 

-1-0.08 

-0.04 

-1-0.05 

0.12 

[Co/Fe] 

I 

-1-0.13 

-1-0.00 

-1-0.01 

-1-0.02 

0.13 

[Ni/Fe] 

I 

-0.03 

-1-0.02 

-1-0.05 

-1-0.01 

0.06 

[Cu/Fe] 

I 

-0.02 

-1-0.03 

-1-0.08 

-1-0.02 

0.09 

[Zn/Fe] 

I 

-0.18 

-1-0.07 

-1-0.04 

-1-0.05 

0.20 

[Sr/Fe] 

II 

-1-0.01 

-0.06 

-1-0.02 

-0.07 

0.09 

[Y/Fe] 

II 

-1-0.09 

-0.01 

-1-0.02 

-1-0.00 

0.09 

[Zr/Fe] 

II 

-1-0.05 

-0.01 

-1-0.02 

-1-0.00 

0.05 

[Mo/Fe] 

I 

-1-0.03 

-1-0.00 

-1-0.05 

-0.01 

0.06 

[Ba/Fe] 

II 

-1-0.16 

-0.02 

-0.08 

-1-0.00 

0.18 

[La/Fe] 

II 

-1-0.10 

-0.02 

-1-0.05 

-1-0.01 

0.11 

[Ce/Fe] 

II 

-1-0.07 

-1-0.00 

-1-0.05 

-1-0.02 

0.09 

[Pr/Fe] 

II 

-1-0.13 

-1-0.02 

-1-0.08 

-1-0.03 

0.16 

[Nd/Fe] 

II 

-1-0.13 

-1-0.00 

-1-0.04 

-1-0.02 

0.14 

[Sm/Fe] 

II 

-1-0.11 

-0.01 

-1-0.03 

-1-0.03 

0.12 

[Eu/Fe] 

II 

-1-0.09 

-1-0.02 

-1-0.07 

-1-0.01 

0.12 

[Gd/Fe] 

II 

-1-0.09 

-1-0.05 

-1-0.04 

-0.01 

0.11 

[Tb/Fe] 

II 

-1-0.10 

-1-0.03 

-1-0.07 

-1-0.05 

0.14 

[Dy/Fe] 

II 

-1-0.07 

-0.02 

-0.02 

-0.02 

0.08 

[Er/Fe] 

II 

-1-0.00 

-1-0.03 

-1-0.00 

-1-0.01 

0.03 

[Tm/Fe] 

II 

-1-0.09 

-0.04 

-1-0.00 

-0.08 

0.13 

[Yb/Fe] 

II 

-1-0.21 

-0.09 

-1-0.04 

-0.08 

0.25 

[Hf/Fe] 

II 

-1-0.12 

-1-0.07 

-1-0.09 

-0.02 

0.17 


Notes— (JTefj = -1-97 K, 51og g = -10.21 dex, Svt = -10.18 km s (5[M/H] = -10.20; 
see Sections 16.II and 16.21 for details. 


are found to overlap the ICarretta et al.l data relatively well 
in the [0/Fe] versus [Na/Fe] plane. 

Our data and that of ICarretta et al.l (l2014l l indepen¬ 
dently exhibit a broad anti-correlation between [0/Fe] and 
[Na/Fe]. This includes stars with [0/Fe] and [Na/Fe] ratios 
consistent with normal halo field stars and stars with de- 
pleted [0/Fe] and enhanced [Na/Fe] ratios. ICarretta et al.l 
calculated that 31 ± 8 per cent of stars in NGC 4833 belong 
to the primordial component (i.e., those with normal [0/Fe] 
and [Na/Fe] ratios). Another 59 ± 11 per cent of stars belong 
to the intermediate component (i.e., those with depleted 
[0/Fe] ratios and enhanced [Na/Fe] ratios), and 10 ± 4 per 
cent of stars belong to the extreme compo nent (i.e., those 
with [0/Na] < —0.9: ICarretta et al.ll200^ ) . Our sample is 
more than three times smaller than that of ICarretta et ahl . 
so we do not attempt to rederive these fractions. The frac¬ 
tion of stars in the primordial component, one-third, is typ¬ 
ical of most globular clusters that have been studied. The 


fraction of stars in the extreme component is larger than 
in most clusters, whic h have typically 2-3 per cent or less 
dCarretta et al.ll2009a^ . 

Figure [8] also reveals correlations between [O/Fe] and 
[Mg/Fe], [Na/Fe] and [Al/Fe], possibly [Na/Fe] and [Si/Fe], 
and possibly [Al/Fe] and [Si/Fe]. Anti-correlations are also 
present among [O/Fe] and [Al/Fe], [Na/Fe] and [Mg/Fe], 
and [Mg/Fe] and [Al/Fe]. These featur es were also present 
in the sample of ICarretta et all ll2014h . Collectively, these 
abundance patterns indicate that proton-capture reactions 
operating at high temperatures occurred in the stars that 
enriched some of the present-day members of NGC 4833. 

ICarretta et al.l (120141 1 identified a bi-modal distribution 
among the Mg and A1 abundances derived from the UVES 
spectra. Our data confirm this separation, which is eas¬ 
ily discerned in the [O/Fe] versus [Al/Fe], [Na/Fe] versus 
[Mg/Fe], [Na/Fe] versus [Al/Fe], and [Mg/Fe] versus [Al/Fe] 
planes in Figure |8| 
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[O/Fe] 

Figure 8. Rel ationships among O, N a, Mg, Al, and Si in NGC 4833. Black symbols mark data from this study, and small gray circles 
mark data from ICarretta et al.l J 2 OI 4 I '). The dotted lines indicate the solar ratios. 


Trends involving [Si/Fe] are more subtle. ICarretta et al.l 
found evidence of such trends in their sample, which can 
be seen by simple inspection of their Figure 8. Similar in¬ 
spection of our Figure |8] only hints that [Si/Fe] may cor¬ 
relate with [Na/Fe] and [Al/Fe] and anti-correlate with 
[Mg/Fe]. We divide our sample into two groups to test 
whether these correlations are significant. The [Mg/Fe] ra¬ 
tios are well-separated into two groups divided at [Mg/Fe] « 
-1-0.3, and we compute the weighted mean [Si/Fe] ra¬ 
tio within each group. In the Mg-rich group of stars, 
{[Si/Fe]) = -1-0.72 ± 0.03; in the Mg-poor group, ([Si/Fe]) = 
+0.76 ± 0.04. The [Si/Fe] ratios are not significantly differ¬ 
ent in the Mg-rich and Mg-poor groups of stars. 

We also compute p-values for the linear correlation coef¬ 
ficient r between [Si/Fe] and other light-element ratios. The 


p-value gives the probability that a random sample of N 
uncorrelated points would yield an experimentally-derived 
value 15 ]r]. For example, the p-value for the [Na/Fe] and 
]Mg/Fe] ratios shown in Figure[8]is 0.0002. The p-values be¬ 
tween [Si/Fe] and [O/Fe], [Na/Fe], ]Mg/Fe], and [Al/Fe] are 
0.645, 0.077, 0.019, and 0.356, respectively. Only the anti¬ 
correlation between [Mg/Fe] and [Si/Fe] is significant at the 
2a level in our data. 


Censored (“survivor”) statistical tests can be applied 
to check whether the [Si/Fe] correlations may be signihcant 
if upper limits are included. We apply the generalized ver¬ 


sion of Kendall’s r rank correlation test to our data using 
the ASURV cod e CRev. l:lLavallev. Isobe. fc Feigelsonlll992l ~l 
as described in llsobe. Feigelson, fc~elsorJ 1 198a ). This test 
evaluates whether a correlation may be present between two 
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14 Ian U. Roederer and Ian B. Thompson 


Table 8. Adopted solar abundances 


Element 

Z 

loge© 

Photospheric (P) 
or meteoritic (M) 

c 

6 

8.43 

p 

N 

7 

7.83 

p 

O 

8 

8.69 

p 

Na 

11 

6.24 

p 

Mg 

12 

7.60 

p 

A1 

13 

6.45 

p 

Si 

14 

7.51 

p 

K 

19 

5.03 

p 

Ca 

20 

6.34 

p 

Sc 

21 

3.15 

p 

Ti 

22 

4.95 

p 

V 

23 

3.93 

p 

Cr 

24 

5.64 

p 

Mn 

25 

5.43 

p 

Fe 

26 

7.50 

p 

Co 

27 

4.99 

p 

Ni 

28 

6.22 

p 

Cu 

29 

4.19 

p 

Zn 

30 

4.56 

p 

Rb 

37 

2.52 

p 

Sr 

38 

2.87 

p 

Y 

39 

2.21 

p 

Zr 

40 

2.58 

p 

Nb 

41 

1.46 

p 

Mo 

42 

1.88 

p 

Ba 

56 

2.18 

p 

La 

57 

1.10 

p 

Ce 

58 

1.58 

p 

Pr 

59 

0.72 

p 

Nd 

60 

1.42 

p 

Sm 

62 

0.96 

p 

Eu 

63 

0.52 

p 

Gd 

64 

1.07 

p 

Tb 

65 

0.30 

p 

Dy 

66 

1.10 

p 

Ho 

67 

0.48 

p 

Er 

68 

0.92 

p 

Tm 

69 

0.10 

p 

Yb 

70 

0.92 

M 

Hf 

72 

0.85 

P 

Ir 

77 

1.38 

P 

Pb 

82 

2.04 

M 

Th 

90 

0.06 

M 


variables when upper limits are found in either variable. We 
find that correlations among [Na/Fe], [Mg/Fe], and [Al/Fe] 
are all highly significant {p < 0.004 in each case). We also 
find that [Si/Fe] does not exhibit a significant correlation 
with either [Na/Fe] or [Al/Fe] if the upper limits are con¬ 
sidered (p > 0.12 in each case). 

We are unable to reproduce correlations between [Si/Fe] 
and any other ratios, except perhaps [Mg/Fe], in our 
data. Our sample siz e is considerably smaller than that of 
ICarretta et al.l (|2014l l. which may explain the difference in 
our results. 


9.3 Potassium through Zinc 

Figure|9]illustrates the relationships among [Mg/Fe], [K/Fe], 
[Ca/Fe], and [Sc/Fe]. Visual inspection of Figure |9] suggests 



j I I^I 

1.0 2.0 
[N/Fe] 


Figure 7. Relationships among C, N, and Na in NGC 4833. 
Black symbols mark data from this study. The dotted lines indi¬ 
cate the solar ratios. 


that [K/Fe] might anti-correlate with [Mg/Fe]. We calculate 
the mean [K/Fe] ratios in the Mg-rich and Mg-poor groups 
of stars. We find ([K/Fe]) = -1-0.36 ± 0.03 for the Mg-rich 
group. We find ([K/Fe]) = -1-0.45 ± 0.03 for the Mg-poor 
group. These ratios are only distinct at the ~ 1.5 ct level, 
which is not significant. The p-value for the linear correlation 
between [Mg/Fe] and [K/Fe] is 0.018, which is significant. If 
[Mg/Fe] and [K/Fe] are anti-correlated, we might also ex¬ 
pect correlations between [K/Fe] and [Na/Fe], [Al/Fe], or 
[Si/Fe]; however, the p-values for the correlations between 
[K/Fe] and other light elements are not significant at the 2a 
level. We conclude that the evidence for an anti-correlation 
between [Mg/Fe] and [K/Fe] is curious but not compelling. 

Hints of an anti-correlation between [Mg/Fe] and 
[Sc/Fe] are seen in Figure [O] but a closer analysis sug¬ 
gests otherwise. Figure [S] demonstrates that the three low¬ 
est [Sc/Fe] ratios are found in three of the four coolest stars 
in our sample, signaling that these low ratios could be an 
artifact of the analysis. This correlation appears weaker if 
only the 11 stars with T^fi > 4500 K are considered. When 
dividing the sample of all 15 stars into Mg-rich and Mg- 
poor groups, we find ([Sc/Fe]) = -1-0.16 ± 0.03 in the Mg- 
rich group and ([Sc/Fe]) = -1-0.25 ± 0.04 in the Mg-poor 
group. If only the 11 stars with TbH > 4500 K are consid¬ 
ered, we hnd ([Sc/Fe]) = -1-0.19 ± 0.01 in the Mg-rich group 
and ([Sc/Fe]) = -1-0.27 ± 0.04 in the Mg-poor group. The 
[Sc/Fe] ratios are not distinct at the 2a level. We also con¬ 
duct a line-by-line differential analysis of the [Sc/Fe] ratios 
in two stars (2-185 and 2-918) with similar stellar param¬ 
eters (Teff = 4670 K and 4696 K, respectively) and differ¬ 
ent [Mg/Fe] ratios ([Mg/Fe] = -1-0.09 and -1-0.40, respec¬ 
tively). An unweighted average of these differentials yields 
5[Sc/Fe] = -1-0.018 ± 0.043 dex, which ind icates no differ¬ 
ence. Finally, we note that ICarretta et al.l (l2014l f found no 
evidence for an anti-correlation between [Mg/Fe] and [Sc/Fe] 
in their data, which are also shown in Figure We conclude 
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[Mg/Fe] 

Figure 9. Relationships among Mg, K, Ca, and Sc in NGC 4833. Black symbols mark data from this study, and small gray circles mark 
data from ICarretta et al.l ll2014l ~). The dotted lines indicate the solar ratios. 


that there is no compelling evidence for variations in the 
[Sc/Fe] ratios within NGC 4833. 

ICohen. Huang, fc Kirbvl(1201 il l. ICohen fc Kirbvl ll2012h . 
and iMucciarelli et al.l 1 I 2 OI 2 I I detected an anti-correlation be¬ 
tween [Mg/Fe] and [K/F e] in the massive, metal-p oor glob¬ 
ular cluster NGC 2419. iMucciarelli et all (l2015l l found a 
smaller, but significant, anti-correlation between [Mg/Fe] 
and [K/Fe] in NGC 2808. These two massive clusters 
are t he only ones where such relations have been identi¬ 
fied (ICarretta et al.l[2013ll . Our results indicate that such 
extreme variat i ons i n K do not occur in NGC 4833. 

I Cohen fc Kirbvl 1 I 2 OI 2 II also detected an anti-correlation be¬ 
tween [Mg/Fe] and [Sc/Fe] in their NGC 2419 data, and 
there is no compelling evidence for a similar relation in 
NGC 4833. 

Figure [To] compares the mean [X/Fe] ratios for elements 
in and near the iron group (loosely considered here as K 
through Zn; 19 ^ ^ 30) in NGC 4833 and normal metal- 

poor halo field stars. The comparison stars are drawn from 
the sa mple of 98 red giants examined by iRoederer et al.l 
(I 2 OI 4 II . Upper limits in the comparison sample have been 
omitted from Figure [10] for clarity. Overall, there is excel¬ 
lent agreement between the mean [X/Fe] ratios for elements 
in the iron group in NGC 4833 and normal halo field stars. 

ICratton, Sneden, fc Carrettal (l2004] i pointed out that 


the [Ca/Fe] in tw o stars in NGC 4833 reported by 
ICratton fc Ortolanil lll989l l is higher by « 0.2-0.3 dex than 
the [Ca/Fe] ratio in other metal-poor globular clusters. The 
ICratton fc Ortolani] r esult w as alread y at odds with t h e data 
of Pilachowski et al. (Il983ll . Neither ICarretta et al.l (l2014li 
nor we reproduce this high [Ca/Fe] ratio in NGC 4833. 


9.4 Neutron-capture elements 


Figure [TT] illustrates the heavy-element abundance distribu¬ 
tions in each star observed in NGC 4833. For comparison, 
Figure [TT] also shows the abundance distribution of the main 
component of the r-process as found in the metal-poor gi¬ 
ant C S 22892-052 (ISneden et al.ll200l . l2009l : IRoederer et al.l 
I 2 OO 9 II . There is generally good agreement among the heavy- 
element abundance distribution from one star to the next in 
NGC 4833. Large variations in the hea vy-element distribu¬ 
tions , like those found in ui Cen (e.g., Norris fc Da Costal 
I 1995 I: ISmith et al.l [ 2 OO 0 IL M2 dYong et al.1 120141) . M22 


I Marino et al.l 20091: 


Roederer et al.l 




or NGC 1851 


I Yong fc Grundah]||2008l :' Carretta et al.ll201 il l are excluded. 


Figure [12] shows the mean heavy element abundance 
distribution for all 15 stars observed in NGC 4833. Three 
templates are shown for comparison. The red line marks the 
main component of the r-process, and this line is identi- 
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Figure 10. [X/Fe] ratios for elements, X, in the iron group. The large black squa re indicates the mean abundances in NGC 4833, and 
the small gray diamonds represent individual red giant halo field stars analyzed by iRoederer et al.l ||2014| ~). The dotted lines indicate the 
solar ratios. 


cal to that found in Figure 1111 The thick gold line marks 
the distribution found in the metal-poor giant HD 122563 
llHonda et al.l l2006l : IRoederer et al.lT2012l i. frequently re¬ 
ferred to as the distribution produced by the weak com¬ 
ponent of the r-process. The long-dashed blue line marks 
the distribution predicte d by the main and strong com¬ 
ponents of the s-pro cess JSneden. Cowan, fc Gallinol l2008l : 
iBisterzo et al.l 1201 il l. These curves illustrate the general 
characteristics of r- and s-process nucleosynthesis and are 
not intended to be rigid representations. For example, 
the abundance distributions that result from s-process 
nucleosynthesis depend on the parameters of stars pass¬ 
ing through the therm ally-pulsing asymptotic giant branch 
phase of evolution (e .g. jBusso. Gallino. fc Wasserburg|l9^ : 
IBisterzo et al.ll2010ll . The mean heavy element abundance 
distribution in NGC 4833 traces the main component of the 
r-process within « 2(j for the elements with Z ^ 57. Ba 
{Z = 56) is enhanced by ~ 0.3 dex relative to the template 
for the main component of the r-process when normalized to 
Eu. Sr, Y, Zr, and Mo {Z = 38, 39, 40, and 42, respectively) 
are enhanced by « 0.9, 0.5, 0.9, and 0.6 dex. 

One possible explanation of the Sr, Y, Zr, Mo, and 
Ba enhancements could be a small contribution from some 
form of s-process nu cleosynthesis. If so, the results of 
IRoederer et al.l (I 2 OI 0 I I suggest that we might expect en- 



Atomic Number 

Figure 12. Average heavy element abundance pattern in 
NGC 4833. Filled squares mark detections, and arrows mark 3cr 
upper limits. The red line, thick gold line, and long-dashed blue 
line represent the main component of the r-process, the weak com¬ 
ponent of the r-process, and the main and strong components of 
the s-process as described in the text. Each of the three curves 
has been renormalized to the mean Eu abundance in NGC 4833. 
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Figure 11. Heavy element abundance pattern for each star observed in NGC 4833. The red lines mark the template of the main 
component of the r-process (CS 22892—052). Filled symbols mark detections, and downward-facing triangles mark 3cr upper limits. 
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hancement of other elements with a substantial s-process 
component in solar system material, such as Yb, Hf, and Pb 
{Z = 70, 72, and 82). No such enhancements are observed. 
Another possible explanation is that the weak component 
of the r-process has contributed to the heavy elements in 
NGC 4833. Halo field stars, like HD 122563, whose heavy 
elements have been produced by this mechanism do not typ¬ 
ically show Ir {Z = 77), which is detected in NGC 4833. Our 
data suggest that the most likely explanation for the heavy 
element abundance distribution in NGC 4833 is a combi¬ 
nation of the main and weak components of the r-process. 
A single r-process with physical characteristics intermedi¬ 
ate between those of the main and weak components of the 
r-process could also be responsible. 

Figure [T3] compares the heavy element abundance dis¬ 
tribution in NGC 4833 with six other metal-poor globular 
clusters. These clusters are selected for comparison because 
large numbers of heavy elements have been studied in multi¬ 
ple red giants within each cluster. The top panel of Figure [T3l 
compares NGC 4833 with three clusters wi th [Fe/H] < 
—2.0: M15 ([Fe/H] = —2.53, R GB stars only; |Sobec k_gt__^ 
[20 i 3), M92 ([Fe/H] = —2.70; iRoederer fc Sned ~ j201llh 
and NGC 2419 ([Fe/H] = —2.06; Cohen et al. 201lll . The 
bottom panel of Figure [T3] compares NGC 4833 with three 
other clusters with ]Fe/H] > — 2.0: M2 ([Fe/H] = —1.68, r- 
only group of st ars: lYong et al, l2014h . M5 ([Fe/H] = —1.40, 
RGB stars only; iLai et al.ll201lil . and M2 2 (]Fe /H] = -1.81, 
r-only group of stars; Roederer et al. 1201 ih . The distri¬ 
butions are normalized to the Eu abundance to eliminate 
differences in the overall amount of heavy elements from 
one cluster to another. Regardless, these differences are 
small: {[Eu/Fe]) = -fO.36 (NGC 4833), -hO.55 (M15), -hO.54 
(M92), -fO.38 (M2), -fO.46 (M5), -fO.35 (M22), and -h0.30 
(NGC 2419). 

A few general trends emerge from this comparison. 
First, the [X/Eu] ratios are consistently sub-solar when “X” 
represents the lighter n-capture elements. The [Y/Eu] and 
[Zr/Eu] ratios are similar in all of the clusters where they 
have been examined. Small enhancements in Sr, Y, Zr, Mo, 
and Ba relative to the main component of the r-process are 
not unusual in metal-poor globular clusters. Second, there 
is a consistent upward trend in the [X/Eu] ratios when “X” 
represents Ba through Sm. Third, there is a consistent down¬ 
ward trend in [X/Eu] when “X” represents Gd through Yb. 
In summary, the heavy element distribution in NGC 4833 
closely resembles that in other well-studied metal-poor glob¬ 
ular clusters. 


9.5 Heavy Element Dispersion in NGC 4833 


In many globular clusters, the heavy elements appear to be 
homogeneous at the limit of observations. M15 is unique 
among all known cl usters in that (1) it has a large range 
(spanning > 0.7 dex. lSneden et al]ll997l : IWorlev et al.ll2013l ~l 
of heavy-element abundances relative to Fe, (2) the star-to- 
star abundance distribution remains relatively constant and 
appears to have originated predominantly via r-process nu¬ 
cleosynthesis, and (3) the heavy-element dispersion shows 
no correlation with the light element disp ersion. These 
charac t eristi cs of M15 have b een v er ified by ISneden et al.l 
lll997l. 200Clh . Preston et al.l (l2006fl. Otsuki et al. 1 2006ll . 
ISobeck et al'l l^2011^ . and Worlev et al] ~( 2013fl . 
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Figure 13. Average heavy element abundance pattern in 
NGC 4833 and six other metal-poor globular clusters. All ratios 
have been normalized to Eu. The three comparison clusters shown 
in the top panel have ([Fe/H]) < —2.0, and the three compari¬ 
son clusters shown in the bottom panel have ([Fe/H]) > —2.0. 
References are given in the text. 


iRoedererl (1201111 reanalyzed literature data and pre¬ 
sented evidence that several other metal-poor clusters (M3, 
M5, M13, and NGC 3201) may also exh ibit a less-extreme 
r-process dispersion. IRoederer fc Sned"^ 1I2OIIII reported an 
r-process dispersion in M92, w hich w a s not found in the 
higher-quality data obtained bv ICohenI (1201 ill . The [La/Fe] 
and ]Eu/Fe] ratios showed a high degree of correlation for 
stars within each of these clusters. This would not be ex¬ 
pected if the La and Eu abundances were independent mea¬ 
sures of a sample with no cosmic dispersion in La or Eu. The 
]La/Eu] ratios were constant, indicating that the distribu¬ 
tions were constant and only the total amount of material 
(relative to, e.g., H or Fe) was changing. In principle, er¬ 
rors in the choice of model atmosphere could also produce 
the observed correlation since the La ll and Eu ll lines form 
similarly. This possibility was dismissed because no corre¬ 
lation was found among other species that might also form 
similarly, including Sc ll and Ti ll. In these clusters, the 
r-process dispersion did not correlate with the [Na/Fe] ratio, 
which was chosen to represent the light-element dispersion 
commonly found in globular clusters. 

We calculate a p-value of 0.001 for the relation between 
[La/Fe] and [Eu/Fe] in 15 stars in NGC 4833. The null hy¬ 
pothesis that [La/Fe] and [Eu/Fe] are not correlated is re¬ 
jected at a > 3a significance level. Figure [T4l illustrates the 
correlations between [Eu/Fe] and eight heavy-element abun- 
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Element "X" 

Figure 15. Slopes in the [X/Fe] versus [Eu/Fe] relationships, 
where “X” represents a given species. Recall that all ratios are 
constructed using ions with ions (e.g., [La Il/Fe ll]) and neutrals 
with neutrals (e.g., [Na l/Fe ij. The vertical dashed line separates 
the light from the heavy elements. The shaded region represents 
the mean ± la of the slopes of the eight heavy elements shown 
to the right of the vertical dashed line. The horizontal dotted line 
marks a slope of zero. 

dance ratios in NGC 4833. These eight pairs of ratios are 
selected because they have been measured in 14 or 15 of the 
15 stars studied. The correlations bear strong resemblance 
to those identified bv iR.oedereil |201^. Of these eight ra¬ 
tios, only [Sr/Fe], [Ba/Fe], and [Nd/Fe] do not correlate with 
[Eu/Fe] at a > 20- level. 

Weighted linear least-squares fits to the data are also 
shown in each panel in Figure 1141 A slope of -|-1 corre¬ 
sponds to a direct correlation, —1 corresponds to a direct 
anti-correlation, and 0 corresponds to no correlation. The 
values of the slopes and la uncertainties are illustrated in 
Figure 1151 The slopes are remarkably consistent and non¬ 
zero, and the mean slope (0.75 ± 0.07, a — 0.19) is illus¬ 
trated by the shaded box in Figure fTSl Only the slope be¬ 
tween [Eu/Fe] and [Sr/Fe] is not strongly positive among 
the elements produced by n-capture reactions. 

Slopes of the relations between [Eu/Fe] and lighter ele¬ 
ment ratios (Z ^ 30) are shown in Figure[T5]for comparison. 
Most of these slopes are consistent with zero, in contrast to 
the slopes between [Eu/Fe] and other elements produced by 
n-capture reactions. There is no coherent, significant pattern 
among them. This suggests that none of the lighter elements 
correlate with Eu or other n-capture elements. 

Figure [16] illustrates [Na/Fe] and eight of the heavy- 
element ratios shown in Figure (T?) No correlations are ap¬ 
parent, and the p-values indicate that none of these correla¬ 
tions (including [Na/Fe] versus [Sm/Fe], not shown) is even 
la significant. The heavy element dispersion is not related 
to the light element dispersion. 

We have shown that the abundances of the n-capture el¬ 
ements in NGC 4833 correlate closely with one another but 
not the lighter elements. Next, we examine whether these 
correlations reflect true abundance dispersions or are an ar¬ 
tifact of our analysis. Atoms in the same state (e.g., ions) 
would respond similarly to an inappropriate model atmo¬ 
sphere, and the resulting abundance distributions would be 
smeared out to lower and higher ratios in a correlated fash¬ 
ion. Table [9] lists the la uncertainty in each [X/Fe] ratio 
that could be attributed to uncertainties in the model at¬ 


mosphere parameters. For example, these values are 0.11 
and 0.12 dex for [La/Fe] and [Eu/Fe]. The standard devi¬ 
ations of the [La/Fe] and [Eu/Fe] ratios are each 0.13 dex. 
The comparable magnitudes of these values suggest that the 
correlated heavy element abundance ratios may not have a 
cosmic origin. 

To demonstrate this point, we perform the following 
test. Using the values presented in Table [9| we adjust the 
Teff values for each star to force the individual [Eu/Fe] ratios 
to equal the mean [Eu/Fe] ratio. We then apply these Tes 
adjustments to the individual [La/Fe] ratios to derive revised 
[La/Fe] ratios for each star. The standard deviation of the 
revised [La/Fe] ratios is 30 per cent smaller than the stan¬ 
dard deviation of the uncorrected ratios. Other rare earth 
elements exhibit similar responses. We conclude from this 
test that random errors in Teu can account for a substantial 
portion of the dispersion, and we infer that random errors 
in the other model atmosphere parameters may have similar 
impact. 

All of the heavy-element ratios shown to correlate in 
Figure [TT] are derived from ionized atoms, so other ele¬ 
ments detected as ions (Sc, Ti, V, and Gr) might exhibit 
similar characteristics. However, Table [9] indicates that the 
[Sc ii/Fe], [Ti ii/Fe], [V n/Fe[, and [Gr ii/Fe] should not be 
expected to behave like, e.g., [La ll/Fe[ or [Eu ll/Fe]. For 
example, (5[Ti ii/Fe]/(5ut = —0.04, while (5[La ii/Fe]/(5ut = 
-1-0.05. Responses to ST^ff and 51og g are also dissimilar. 
Even though both ratios are derived from lines of ionized 
atoms, the line strengths and excitation potentials are dif¬ 
ferent. The excitation potentials of Sc ll, Ti ll, V ll, and 
Gr II lines used range from 1.08 to 4.07 eV with a median 
of 1.50 eV, while the excitation potentials of the heavy rare 
earth elements’ lines range from 0.00 to 1.38 eV with a me¬ 
dian of 0.32 eV. This difference is significant. The [Sc ii/Fe], 
[Ti ll/Fe[, [V ii/Fe[, and [Cr ll/Fe] ratios are not a good 
control group for [La ll/Fe] or [Eu ll/Fe], as was assumed by 
iRoederej 1I2OIIII . 

We conclude that there is no compelling evidence of a 
cosmic origin for the correlations among pairs of n-capture- 
element ratios in NGC 4833. The detected correlations likely 
result from random errors in the model atmosphere parame¬ 
ters. We sugg est that the corr elations identified among liter¬ 
ature data bv iRoederer] ll201lji may also be the result of simi¬ 
lar random errors in the model atmosphere parameters. The 
r-process dispersion in M15, however, is too large (> 0.7 dex) 
to b e explained entirely by this ph enomenon, and we note 
that iTsuiimoto fc Shigevamal (l2014l ~) have recently offered a 
theoretical explanation of its origin. 

9.6 The nuclear chronometer in NGC 4833 

We derive upper limits on the Th abundance based on the 
non-detection of the Th ll line at 4019 A. This line origi¬ 
nates from radioactive ^®^Th, which can only be produced 
via r-process nucleosynthesis. This upper limit sets a lower 
limit on the age of the r-process material in NGC 4833. We 
compare the Th abundance to Eu, a well-measured, stable 
element presumably produced by the same nucleosynthesis 
channel. The lowest Th/Eu ratio found in NGC 4833, in star 
2-1664, is loge(Th/Eu) < —0.47 ± 0.09. If we assume the 
T h and Eu were produc ed in the ratio reported in Table 9 
of iRoederer et all (l2009l L this implies an age greater than 
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Figure 14. Correlations among heavy element ratios in NGC 4833 for elements measured in 14 or more stars. The red lines mark 
the weighted least-squares fits. Dotted lines mark the solar values. We have checked, and confirmed, the weighted fit for [Ba/Fe] versus 
[Eu/Fe], which appears to deviate from what might be expected from a by-eye fit. The uncertainty on this fit is considerably larger than 
for the other rare earth elements, as can be seen in Figure fTSl 
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Figure 16. Relationsh ips among heavy elem ents and [Na/Fe] in NGC 4833. Black points mark our derived abundance ratios, and gray 
circles mark those from ICarretta et al.l ||2014|') . Dotted lines mark the solar values. No correlations are apparent. 


6.4 Gyr. This value is not very constraining, but it implies 
that the r-process material in NGC 4833 is not young. 

Alternatively, we could assume that NGC 4 833 is old 
([Melbourne et ai1l2000l : iMarm-Franch et al.ll200^ 'l and treat 
the initial Th/Eu production ratio as a free parameter. An 
assumed age of 12 ± 1.5 Gyr implies an initial production 
ratio of loge(Th/Eu) = —0.21 ± 0.09. This implies that 
the r-process material was not produced in an “actinide 


boost” event (cf. iHill et ^l2002l : ISchatz et al.ll2002h . which 
would require an initial production ratio of log£(Th/Eu) = 
-1-0.11 ± 0.07 to yield an age consistent with the prototype 
actinide boost halo star CS 31082-001. 

To the best of our knowledge, no actinide boost has 
yet been observed in any globular cluste r star. ^^^Th ha s 
been detected in globul ar clusters M5 dLai et al.l 1201 il l. 
M15 (ISneden et ahll^OOll . M22 (Roederer et al. l201lh . and 
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Figure 17. Comparison of globular cluster ages derived from 
isochrone fitt ing and radioactive decay o f The relative ages 

presented bv lMarm -Franc h et al.l ll2009h are based on stellar evo¬ 
lution models bv lDotte^^aP ll2007h and have been normalized 
to an absolute age of 13.0 Gyr. The dotted line marks the 1:1 
correspondence. 

M92 (| Johnson fc Bolt^ |2QQ1^ ■ We also report an unpub¬ 
lished upper limit on in M2, loge(Th/Eu) < —0.42, 

based on spectra obtained by I.U.R. and published in 
lYong et al.l hoidll . Ages calculated from the loge(Th/Eu) 
ratios are consistent with ages deduced by isochrone htting 
dMan'n-Franch et al.ll200^ ~). as shown in Figure [170 The un¬ 
certainties are substantial, but we would expect to End val¬ 
ues in the lower right corner of Figure[T7]if an actinide boost 
is present in any of these six clusters at the level found in 
CS 31082-001. No detections are found here. 

The astrophysical site(s) responsible for r-process nu¬ 
cleosynthesis are still debated. Non-detection of the actinide 
boost in globular clusters may provide a new environmental 
test of r-process nucleosynthesis models. The sample of clus¬ 
ters where ^^^Th has been studied is still small, however, so 
new data are of great importance. 


10 CONCLUSIONS 

We have collected new high-resolution spectroscopic ob¬ 
servations of 15 giants in the metal-poor globular cluster 
NGC 4833. We derive stellar parameters for these stars and 
examine up to 354 lines of 48 species of 44 elements in each 
star. We detect 43 species of 39 elements and present upper 
limits derived from the non-detections of the others. Over¬ 
all, the composition of stars in NGC 4833 appears relatively 
normal for a metal-poor globular cluster. 

We derive a mean metallicity of ([Fe/H]) = 

^ lYong et al.l ll2008al f9i detected Th in M4 and M5, but the high 
[Th/Fe] ratio in M5 derived by lYong et ahl was not confirmed by 
iLai et J ll201lh . We refrain fr om drawing any conclusions from 
the [Th/Fe] ratios presented by lYong et all 


—2.25 ± 0.02 {a = 0.09) from Fe l lines and ([Fe/H]) = 
—2.19 ± 0.013 (cr = 0.05) from Fe ll lines. These uncer¬ 
tainties support the results of ICarretta et al.l |[20l3) that 
there is no internal Fe dispersion within NGC 4833. Our de¬ 
rived mea n metallicity is 0.20 ± 0.02 dex lower than that 
derived bv ICarretta et aP for four stars in common. We de¬ 
rive [Fe/H] = —2.19 ± 0.04 for NGC 4833 on the differen- 
tial gl o bular cluster scale established bv iKoch fc McWilliamI 
(l2008l . I2OIIII relative to the K-giant Arcturus, which sup¬ 
ports our lower metallicity. 

Our data reveal the abundance variations among 
[O/Fe], [Na/Fe], [Mg/Fe], [Al/Fe], and possibly [Si/Fe] com¬ 
monly found in globular clusters. There are bi-modal dis¬ 
tributions among the [Na/Fe], [M g/Fe], and [Al/Fe] ra - 
tios. Our results reaffirm those of ICarretta et ahl (l2014ll . 
We also reproduce a hint of an anti-correlation between 
[Si/Fe] and [Mg/Fe], but our data do not reveal statistically- 
significant correlations between [Si/Fe] and either [Na/Fe] or 
[Al/Fe]. We su spect the discrepancy between our data and 
ICarretta et al.l is a result of the smaller sample size exam¬ 
ined by us. 

Quantitative measures of potential anti-correlations be¬ 
tween [K/Fe] and [Mg/Fe] and between [Sc/Fe] and [Mg/Fe] 
reveal that neither is significant. Furthermore, neither 
[K/Fe] nor [Sc/Fe] correlate with [Na/Fe], [Al/Fe], or [Si/Fe] 
at the 2 ct level. Our data demonstrate that NGC 4833 does 
not possess the extreme relations among [Mg/Fe], [K/Fe], 
and [Sc/Fe] found in the globular clusters NGC 2419 and 
NGC 2808. We find no dispersion among any of the other 
iron-group elements in NGC 4833. 

We detect up to 20 n-capture elements in NGC 4833, 
and we place upper limits on four others. Abundances of 
the rare earth elements are consistent with r-process nucle¬ 
osynthesis. Sr, Y, Zr, Mo, and Ba are enhanced relative to 
the scaled solar r-process distribution, but low levels of Yb, 
Hf, and Pb suggest that this is not due to enrichment by 
the s-process. Instead, the weak and main components of 
the r-process may be responsible. The heavy-element abun¬ 
dance distribution in NGC 4833 closely resembles that found 
in M5, M15, M92, NGC 2419, and the r-process-only stel¬ 
lar groups in M2 and M22. There is no correlation between 
the n-capture elements and the light-element variations in 
NGC 483 3, which reaffirms an d expands upon results ob¬ 
tained bv ICarretta et ah! ll2014ll . Our analysis of correlations 
between [La/Fe], [Eu/Fe], and other n-capture elements in 
NGC 4833 reveals that they are likely due to random er¬ 
rors in the stellar parameters, and we conclude that there 
is no compelling evidence for cosmic dispersion among the 
heavy elements in NGC 4833. This cluster does not possess 
wide star-to-star variations like those observed in to Gen, 
M2, M15, M22, or NGC 1851. 

Upper limits derived from the non-detection of the 
Th II line at 4019 A indicate that the r-process material 
in NGC 4833 was not formed in an actinide boost event. 
No other clusters that have been studied show evidence of 
an actinide boost, either. The potential non-detection of the 
actinide boost in globular clusters presents a new environ¬ 
mental test of r-process nucleosynthesis models that should 
be investigated further. 
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